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We present here a phenomenological cosmological model under perfect fluid distribution with a 
stifle equation of state p = p- The erstwhile cosmological constant is assumed to be a time dependent 
variable, i.e., A = A(t) in our study. It has been shown that the estimates of different cosmological 
parameters from this model are in good agreement with the experimental results, especially 13.79 
Gyr as the age of the universe is quite satisfactory. The behavior and relation of A-stiff fluid model 
with dust, viscous fluid and variable G have also been investigated in detail. 



PACS numbers: 98.80.Jk, 98.80.Cq 

"The problem of the age of the observable uni- 
verse - the time back to the 'big bang' - has puzzled 
laymen, churchmen and scientists for many years." 
W. A. Fowler (in The Age oj the Observable Universe) 0| 

1. Introduction 

The cosmological picture that emerges after the obser- 
vations of SNe la by HZT and SCP teams reveals that at 
present we are residing in an accelerating universe |^ 
whose geometry is Euclidean in nature Q . This speeding 
up of the universe started about 7 Gyr ago and some 
kind of repelling force, termed as dark energy, is sup- 
posed to be responsible for catapulting the once decel- 
erating Universe into an accelerating one. Using various 
observational data it has also been possible to pin down, 
to some extent, the ranges of various parameters of the 
Universe at the present era. For instance, the acceptable 
limit of Ho, the Hubble parameter at present, is 72 ± 8 
Km s~^Mpc~^ 1^ and while taking into consideration 
of various ranges as proposed by different workers from 
varying standpoints the currently accepted observa- 
tional value of the age of the universe becomes 14 ± 0.5 
Gyr ^ (an elaborated data sheet has been produced in 

Many variants of A, a dark energy representative, have 
been proposed by different workers to account the cosmo- 
logical consequences. Since, supernovae la observations 
predict a small value of A at the present epoch , now- 
a-days A with dynamical character is preferred over a 
constant A This has been assumed to vary with 

time so that decrease of A from a large initial value to 
its present small value can easily be realized. 

Mainly, in the theoretical point of view, the possibility 
of a nonzero and varying A came into picture in connec- 
tion to the age problem of the universe. This is because 
of the fact that for large A the age of the universe can, in 



principle, become infinite [T2|. Even with the time- 
decaying A, it is seen that the present age of the uni- 
verse either suffers from low-age problem or it is as 
large as 27.4 ± 5.6 Gyr In connection to the low 

age we can mention the results of Overduin and Cooper- 
stock 
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where for the Hubble parameter Hq = 73 ± 10 
kms~^Mpc^^ the age of the universe in its lower limit is 
5.4 Gyr! The situation improves slightly for lower value 
of i?o(= 55 ± 10 kms"iMpc"i) (l|l when it becomes 7.6 
Gyr only. This values are much lower than the estimated 
globular cluster ages which are thought to be in the range 
9.6 Gyr g to 13±3 Gyr [|l| and very recent observations 
claim that it is most probably 12.5 ± 1.2 Gyr (l|, 0. 
However, in the standard model of a flat A-dominated 
universe the age of the universe is 13.7 ± 0.2 Gyr as ob- 
tained by Spergel et al. |l^ whereas Kunz et al. [|l9j in the 
A-CDM case find out value as 13.55 ±0.26. For a slightly 
closed A-CDM universe with = 6616:1 kms'^Mpc"^ 
Tegmark et al. estimate an age of 14.llJ'g Gyr. 

Under this background, in the present article we have 
investigated three types of widely used ^ phe- 
nomenological forms of kinematical A, viz., A ^ (d/a)^, 
A ^ [a /a) and A ^ p, where a is the scale factor of the 
Robertson- Walker metric and p is the energy density of 
the universe (a detail account of these phenomenological 
A models is available in ref. [^). 

As a result of the present investigations under the dif- 
ferent phenomenological models as mentioned above we 
have found out the equivalence of the A-models through 
the parameters a, (3 and 7 involving in the models. Un- 
der the perfect fluid distribution, specially with a stiff 
equation of state p = p, the age and the other physical 
parameters have been estimated. Surprisingly, the age 
13.79 Gyr is in good agreement with the present avail- 
able experimental results. The behavior and relation of 
stiff fluid to dust and viscous fluid also have been criti- 
cally investigated. It is also observed that Variable A and 
variable G models show some unique common features in 
connection to the present accelerating universe. 

The article has been organized as follows: in Sec- 
tion 2 the standard general relativistic Einstein fleld 
equations are presented with the inclusion of time 
varying cosmological constant, viz., A = A(t), in the 
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energy-momentum tensors. Section 3 deals with the 
solutions of various phenomenological models within 
this framework of modified general relativity. In the 
Section 4 some salient features of the present model are 
discussed. Section 5 presents concluding remarks on the 
status of the stiff fluid and also the problems related to 
the age of the universe. 

2. Einstein Field Equations 

Let us consider the Robertson- Walker metric which is 
given by 



-aity 



dr^ 



1 — kr^ 



r^{de^ 



(1) 



where k is the curvature constant which takes the specific 
values —1,0 and -1-1 respectively for open, flat and close 
models of the universe and a is the cosmic scale factor as 
mentioned earlier. 

With the assumption that the so-called cosmological 
constant is time dependent here, viz., A — A(i) and c, 
the velocity of light in vacuum is assumed to be unity in 
relativistic units, the Einstein field equations 
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for the above spherically symmetric metric (1) yield, re- 
spectively, the Friedmann and Raychawdhuri equations 
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where G, p and p are, respectively, the gravitational con- 
stant, matter-energy density and fiuid pressure. 

Also, the energy conservation law can be written as 



, .a SttG . A 

inG{p + p)- = —p - -. 

a 6 A 



(5) 



Let us now choose the barotropic equation of state in the 
form 



P = wp, 



(6) 



where w is a constant, known as the equation of state pa- 
rameter and can take the values 0, 1/3 and 1 respectively 
for the pressureless dust, electromagnetic radiation and 
stiff or Zel'dovich fluid. 

Using this barotropic equation of state (6) , in a straight 
forward way, the equation (4) can be written as 
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a 6 6 



(7) 



After performing differentiation of equation (3) with re- 
spect to time and using equations (4) - (7) for elimination 
of p, one obtains the following equation 
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This is the general equation for our investigation in the 
presence of the curvature parameter k. 

3. Solutions 

At this point let us consider the ansatz A — 3Q;(a/a)^, 
where a is a constant. Now, inflation theory predicts 
and CMB detectors such as BOOMERanG @, H, ||, 
MAXIMA 11,11,11, DASI ||, CBI [|o) and WMAP 
pl| , |3^ confirm that the universe is spatially flat. There- 
fore, for the flat universe where k = Q, the equation (8) 
reduces to 
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0. 



(9) 



The general solutions of the above equation (9) can be 
given as 

a(i)-C'ii2/3(i-")(i+-), 
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where Ci is an integration constant and < a < 1 for 
physical validity. 

Therefore, for stiff fiuid distribution (w = 1) the pro- 
portional relation of the cosmic scale factor to the cos- 
mic age becomes in a similar way, it can be 
shown that for the ansatz A — [3{a/a) and A — SnG^p 
the cosmic scale factor a{t) is, respectively, proportional 
to i(/3-i)/(/3-3) and t('>'+i)/3, where /? and 7 are another 
two constants like a. For all these three models, as are 
evident from equations (11) and (12), the cosmic matter- 
energy density p(t) and the cosmological parameter A{t) 
both follow an inverse square law with t. 

Now, in absence of any curvature, the cosmic matter- 
energy density parameter flj^i— SirGp/SH^) and cosmic 
vacuum-energy density parameter flA{= A/Si/^) are re- 
lated by 



1. 



(13) 



This result is consistent with the current constraint on 
cosmic density parameters |]3^ . Expressing the solutions 
of (8) for stiff fluid in terms of fim, it has been found that 
for all the three A-models, a{t) is proportional to t^^^^"^ 
whereas the expressions for p{t) and A{t) are identical. 
Moreover, a, (3 and 7 are related to Vim and JIa by 
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Equation (14), with the help of equation (13), enables us 
to have interrelation between a, (3 and 7 as follows 



2(3 



7 



3(/3-l) 1+7^ 



(15) 



This clearly shows that the three forms, viz., 
A — 3a(a/a)^, A — (3{d/a) and A = SnG^p, are 
equivalent and the three parameters a, (3 and 7 are 
interconnected by the relation (15). Therefore, it is 
possible to search for the identical physical features of 
others if any one of the phenomenological A model is 
known. 

4. Physical Features 

4.1. Age of the universe from the stiff fluid 
model 

The most striking and puzzling feature of this stiff fluid 
model is the fine agreement of the physical parameters, 
especially the present age of the universe, with the obser- 
vational results. To get an explicit expression for age let 
us differentiate equation (10) and divide it by the cosmic 
scale factor a, which ultimately yields 



t 



•i{\-a){l+w)H' 



(16) 



In the case of stiff fluid [w — 1) the above equation (16), 
by the use of the equations (13) and (14), reduces to 



t = 



1 



(17) 



Recent measurements indicate that the value of rimOi 
the cosmic matter-energy density parameter at present, 
is = 0.33 ± 0.035 @. Assuming n^no = 0.295, 

the lowest acceptable limit of this range, and Hq = 
72kms~^Mpc~-^, the present age of the universe, to{= 
l/357„iO-ffo) is found to be 13.79 Gyr which is in excel- 
lent agreement with the age of the universe as calculated 
on the basis of WMAP data [|| as well as the current 
accepted age mentioned earlier. A comparison between 
the ages for other values of lo shows a situation which 
goes in favor of stiff fluid (TABLE I). 

Moreover, the values of po{— l/247rG'ri„ioio^) and qo, 
respectively the present values of the cosmic matter- 
energy density p and the deceleration parameter 5, 
as calculated on the basis of 0^0 = 0.33 ± 0.035 and 
Hq ~ 72kms^^Mpc^^, confirm that the present models 
support the idea of an open as well as accelerating uni- 
verse [|4|. Also, the value of Ao(= (1 — Qmo) / '^^mo'^to'^) 
^ 10~^~^ is in agreement with the current status of A 
with a small value psj . 

4.2. The enigma of the stiff fluid model 

It is worthwhile to mention that from the solutions of 
equation ( 9), one may arrive at the solutions obtained 
by Arbab ]3q| for the dust {cu = 0) case under the ansatz 
A = 3Q;(d/ap. This is possible via the transformation 



relation in between a and (3 in the form a = f3/3{(3 — 2). 
By the use of this one can easily arrive at the equations 
(10) - (12) which are the same as those of his equa- 
tions (8) - (10). However, in the case of stiff fluid, the 
expected transformation relation can be obtained from 
our equation (15) and his equation (18) which yields 
the value for (3 = 3. It can be seen that this immedi- 
ately makes the whole set of Arbab's equation to blow 
up. This, therefore, suggests that such transformation 
is forbidden which is also true in his own case as this 
makes the solution set singular. However, exploiting 
this transformation relation a = (3/3{(3 — 2) in the so- 
lution of Majernfk |3^ it has been shown by Arbab |3^] 
that a can be expressed as having the form k/(1 + k), 
where k = /3/2(/3 — 3) is a free parameter. This yields 
a = -wq and wq = Pq/pq (-1 < wq < 0), being the 
Quintessence equation of state, Arbab argued that 
Wq is nothing but the vacuum energy parameter. One 
can easily find out that the k of Majernfk ||3^ is nothing 
but our 7 under the ansatz A = SirG^p and follows the 
same relation as shown in equation (15). It suggests that 
stiff fluid has a deep connection to dark energy. 

It is also interesting to note that by assuming n = Q\ 
(this assumption is not unreasonable since, 2/3 < n < 1) 
in the equation (27) of the bulk viscous model of Arbab 
psf with variable G, we get the same expression for the 
age of the universe, viz., to(= ^/^flmoHo) in our model. 
Therefore, it is possible to obtain the same value of to for 
the dust case of Arbab as that obtained in the present 
stiff fluid model. Now, it is known that bulk viscosity 
is associated with the inflationary universe scenario and 
after the decoupling of neutrino in the early universe mat- 
ter did behave like viscous fluid |^ . This bulk viscosity, 
therefore, is similar to a variable cosmological term with a 
repulsive pressure. In the present investigation, variable 
A-based stiff fluid model with constant G and Arbab's 
viscous model with variable A and variable G produce 
the same result so far as the present age of the universe 
is concerned. This again suggests about the underlying 
relation between the stiff fluid and the viscous fluid which 
has a feature of dark energy. 

In this context we would like to add here that 
Chakrabarty and Pradhan |4C[| , with a time-dependent 
gravitational constant G, have shown that the cosmolog- 
ical constant A varies as as in the present stiff fluid 
case with constant G. So, there seems to exist some 
connection between the models with variable G and the 
present stiff fluid-filled Zel'dovich universe. 

5. Concluding Remarks 

Although in the present cosmological scenario it is cus- 
tomary to choose the value of the barotropic index lu as 
zero (for pressureless dust) and sometimes as 1/3 (for ra- 
diation), yet stiff fluid model which refers to a Zel'dovich 
universe have been selected by some authors for vari- 
ous situations such as cold baryonic universe ]4l]] , early 
hadron era Ell , scalar field fluid lEsIl , for the relativistic 
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TABLE I: Age of the universe from different uj models 
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12.43 
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situation prevailing during the early stages of the uni- 
verse and LRS Bianchi I cosmological models . There 
are recent applications and claims for equation of states 
in the various astrophysical realm, e. g. in neutron star 
RX J1856-3754 and hyperon stars |Q which are very 
close to the stiff fluid limit (for some more astrophysical 
as well as cosmological applications, see and refs. 
therein) . In this connection we are interested to mention 
the work of Buchert |^ where it has been shown that 
in the spatially averaged inhomogeneous cosmologies the 
averaged equations show that the averaged scalar curva- 
ture must generically change in the course of structure 
formation and that an averaged inhomogeneous perfect 
fluid, in some cases, act like a free scalar field source 
which can be modelled by a stiff fluid. 

It is seen that the present age of the universe as 
calculated in the standard Friedmann model as well 
as in some other models either suffers from low-age 
problem or it is as large as 27.4 ± 5.6 Gyr in the 
'favored' dust case. In that respect the present stiff fluid 
model with variable cosmological constant demands 
some attention because the age of the model fits well 
with the modern accepted limit. But, the result of this 
A-stiff fluid model is very much surprising in the sense 
that we do not know how this value falls within the 
range of the present values of the age of the universe as 
obtained in the different A-CDM models. It is already 
discussed in the earlier subsection that there exists some 
kind of underlying relation in between either stiff fluid 



and variable G model or stiff fluid and viscous fluid 
which has a feature of dark energy or with the both. If 
there is some relation between them then what is the 
possible mechanism - that is also not understood. One 
possibility may be that the effect of the stijfness of the 
stiff fluid through the process of evolution via radiation 
to the present epoch of matter-dominated universe is 
too weak to perceive and hence dark energy in terms 
of variable cosmological constant does play its definite 
role for the A-stiff fluid model. Thus, here also we are 
facing a kind of age problem of the universe which is 
really intriguing! In this respect the comment made by 
Born, on different context, seems appropriate to quote: 
"Whether one or the other of these methods will lead 
to the anticipated "world law" must be left to future 
research" 
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